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Abstract 

We consider the partial difference equations of the Adler-Bobenko-Suris classification, 
which are characterized as multidimensionally consistent. The latter property leads naturally 
to the construction of auto-Backlund transformations and Lax pairs for all the equations 
in this class. Their symmetry analysis is presented and infinite hierarchies of generalized 
symmetries are explicitly constructed. 

1 Introduction 

It is well known that, integrable differential equations, like the Korteweg-de Vries (KdV), the sine- 
Gordon and the nonhnear Schrodinger, have many properties in common. They can be written 
as the compatibility condition of a Lax pair, which plays a crucial role in solving the initial value 
problem by the inverse scattering transform. They admit auto-Backlund transformations, which 
allow us to construct new solutions from known ones [Ij. 

Other properties arise from the symmetries of these equations. Specifically, they admit infinite 
hierarchies of generalized symmetries and, consequently, infinite conservation laws [23j. Also, they 
reduce to Painleve equations and result from reductions of the Yang-Mills equations [T5^ . All the 
above mentioned characteristics may be considered as criteria establishing the integrability of a 
differential equation. 

Analogous characteristics seem to be in common among integrable difference equations defined 
on an elementary square of the lattice. In the most well known cases, these equations are char- 
acterized by their "multidimensional consistency". For a two dimensional equation, this means 
that, the equation may be imposed in a consistent way on a three dimensional lattice, and, con- 
sequently, on a multidimensional one. This property incorporates some of the above mentioned 
integrability aspects. Specifically, the consistency property provides the means to derive algo- 
rithmically Backlund transformation and Lax pair of the difference equation under consideration 

mm- 

Adler, Bobenko and Suris (ABS) used multidimensional consistency as the key property char- 
acterizing integrable difference equations to classify integrable scalar difference equations [3J S] ■ 
The equations emerged from the classification split into the H and Q list and comprise, apart 
from known ones, a number of new cases. The subsequent study of the resulting equations has 
led to construction of exact solutions [7l [8] , Backlund transformations '6' , symmetries [24l [30] 
and conservation laws [25] . 
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The ABS equations wiU be used as the illustrative example to investigate the similarities 
between discrete and continuous integrable equations, as described above. Specifically, we will 
show how the auto-Backlund is inherited by the multidimensional consistency to these equations. 
Using the well known relation between such transformations and Lax pairs, |10j , we will derive 
Lax pairs for the ABS equations. 

The second direction of our investigation will be the symmetries of the ABS equations. In 
particular, we will prove that they admit infinite hierarchies of generalized symmetries, which are 
constructed inductively using linear differential operators. The latter can be regarded as recursion 
operators for the higher members of the ABS class, i.e. equations Q3 and Q4. 

The paper is organized as follows. In Section 2 we introduce the notation and present the 
background material on symmetries of difference equations. The next section is devoted to a 
short description of two classes of difference equations, one of which contains the ABS equations. 

Section 4 deals with the integrability aspects of the equations under consideration, and, in 
particular, with their auto-Backlund transformations and Lax pairs. Section 5 describes the 
symmetries of the ABS equations and how they can be used effectively to construct solutions, as 
well as infinite hierarchies of generalized symmetries. 

2 Notation and preliminaries on symmetries of difference 
equations 

We first introduce the notation that it will be used in what follows. Also, we present those 
definitions on symmetries of difference equations that will be used in the next sections. For details 
on the subject, we refer the reader to the clear and extended review by Levi and Winternitz [14j . 

A partial difference equation is a functional relation among the values of a function m : Z x Z — > 
C at different points of the lattice, which may involve the independent variables n, m and the 
lattice spacings a, as well, i.e. a relation of the form 

^(n,m,U(o,o),U(i,o),'"(o,i), ■ • • = 0. (1) 

In this relation, is the value of the function u at the lattice point (n + i, rn + j), e.g. 

■"(0,0) = u{n,m), U(i^o) = u{n+l,m), U(o4) = u{n,m+l), 

and this is the notation that we will adopt for the values of the function u from now on. 

The analysis of such equations is facilitated by the introduction of two translation operators 
acting on functions on Z^, defined by 

(^i''^^^) = W(fe,o) , (^if ^u) = "(o,fe) , where fc e Z , 

respectively. 

Let G be a connected one-parameter group of transformations acting on the dependent variable 
W(o,o) of the lattice equation ([T]) as follows 

G : M(o,o) — ' "(0,0) = '^{n,m,U(o^oy,e) , e e R. 

Then, the prolongation of the group action of G on the shifted values of u is defined by 

G'*"^ : ("(i,j)) — > = $(n + i,m + j,U(jj);e)) . (2) 
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The transformation group G is a Lie point symmetry of the lattice equation ([T]) if it transforms 
any solution of ([1]) to another solution of the same equation. Equivalently, G is a symmetry of 
equation ([1]), if the latter is not affected by transformation i.e. 

m, U(o,o),U(i^o),'«(o,i)> ■ • • = 0. 

In essence, the action of the transformation group G is expressed by its infinitesimal generator ^ 
i.e. the vector field 

X = i?(n,m,U(o,o))9„,o_o) , 

where R(n, m, U(o,o)) is defined by 

i?(n,m, U(o,o)) = ^$(n,m,ii(o,o);e) 

e=0 

and is referred to as the characteristic. 

Using the latter, one defines the fc-th order forward prolongation of x, namely the vector field 

k k — i 
i=0 j=0 

This allows one to give an infinitesimal form of the criterion for G to be a symmetry of equation 
([T]). It consists of the condition 

x^*") {S" {n, m, U(o,o), U(i,o), "(04)' ■ • ■ ; a, /?)) = . (3) 

which should hold for every solution of equation ([1]) and, therefore, the latter should be taken 
into account, when condition ([3]) is tested out explicitly. 

Equation ^ delivers the most general infinitesimal Lie point symmetry of equation ([1]). The 
solutions of ([3]) determine the Lie algebra g of the corresponding symmetry group G and the 
latter can be constructed by exponentiation: 

$(n,TO,U(o,o);e) = exp(ex) U(o,o) ■ 

On the other hand, one may consider a group of transformations T acting, not only on the 
dependent variable u, but on the lattice parameters a, (3 as well. This leads to the notion of the 
extended Lie point symmetry. The infinitesimal generator of the group action of F is a vector 
field of the form 

V = R{n,m,U(^o,o))9u^o,o) + (,{n,m, a, (3) da + C,{n,m,a, (3) dp , 

and the infinitesimal criterion for a connected group T to be an extended Lie point symmetry of 
equation ^ is 

v^*") [S (ri,m,U(o,o),U(i,o),W(o,i), • ■ • ;a,/3)) = 0, 

where 

k k—i 

^^'^ = E E {^n^ ° "(0^0)) ^"(..) 

+ £_{n,m,a, (3) da + C{n,m,a, (3) d/3 
is the fc-th forward prolongation of the symmetry generator v. 
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By extending the geometric transformations to the more general ones, which depend, not 
only on n, m and M(o,o)j but also on the shifted values of u, we arrive naturally at the notion of 
the generalized symmetry. For example, a five point generalized symmetry may be given by the 
vector field 

V = i?(n,m,M(o,o),'"(i,o)7'«(o,i),W(-i,o):'"(o.-i))5«(o_(„ , 
while an extended five point generalized symmetry is given by the vector field 

V = i?(n,TO,'U(o,o),W(l,0):"(0,l),W(-l,0):"(0 -1))5«(0.0) 

+ £_{n, m, a, (3) da + ({n, m, a, [3) dp , 

respectively. 



3 Classes of difference equations 

In this section wc introduce two classes of lattice equations, which will be denoted by ^ and , 
respectively. These families contain equations involving the values of a function u at the vertices 
of an elementary quadrilateral, as shown in Figure [T] 
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Specifically, the equations belonging in these classes are autonomous and have the form 

Q(w(o,o),'«(i,o),'"(o,i),U(i,i);a,/3) = 0, (4) 

where the function Q satisfies two basic requirements: it is affine linear and depends explicitly 
on the four indicated values of u, i.e. 

9n(i_j.)Q(M(0,0),W(l,0),'"(0,l),U(l,l);a,^) 7^ 

and 

^«(i,,)Q("(0,0):'"(l,0):"(0,l):'«(l,l);Q;,/3) = 0, 

where i, j — 0, 1. 

Further restrictions on the function Q provide the conditions which separate the members of 
the above classes. Specifically, the members of ^ possess the D4-symmetry, i.e. the function Q 
satisfies the following requirement: 

^("(o.o), "(1,0), "(0,1): ""(1,1); a, /3) = eQ(u(o,o), "(o,i), '"(1,0), "(i,i); /3, «) 

= cr<3("(i,o),"(o.o):'«(i.i),M(o.i);a,/3) , 

(5) 
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where e = ±1 and a = ±1. 

On the other hand, equations belonging in possess the symmetries: 

Qiu{o,o),U(ifi),U(^a,i},U(iX)) = T(3(u(i,o),"(o.o),W(i,i),'"(o,i)) 

= r'Q(u(o,i),U(i4),U(o,o):"(i,o)) , (6) 

where r = ± 1 and r' = ±10. 

It is easy to prove that, actually contains To this end, it is sufficient to prove that, 
every equation in the class ^ also satisfies requirement (|6|). Indeed, supposing Q{u, x,y, z; a, P) 
satisfies ([5]), then using the latter relations one easily finds that 

Q{u,x,y,z;a,(3) = a Q{x,u, z,y; a, (3) = a Q{y, z,u,x;a, (3) , 

i.e. the function Q also satisfies ^ with t — t' — a. 

The affine linearity of Q implies that one can define six different polynomials in terms of Q 
and its derivatives [21 HI [30] , four of them assigned to the edges and the rest to the diagonals of 
the elementary quadrilateral where the equation is defined, as illustrated in Figure [2] 

Specifically, each polynomial hij is defined by 

hij — hj I . — Q Q ,ij Q ,i Q .j ^ ^ ~/~ j ^ ^7 j — 1; ■ ■ ■ 1 '^7 

where Q^i denotes the derivative of Q with respect to its i-th argument and Qjj the second order 
derivative Q with respect to its i-th and j-th argument. It is a bi-quadratic polynomial in the 
values of u assigned to the end-points of the edge or diagonal at which it corresponds. Finally, 
the relations 

/112/134 = /113/124 = /114/123 (7) 

hold in view of the equation Q = 0. 

The imposed symmetries of Q yield further restrictions on hij. In particular, taking into 
account the symmetries on Q for each class separately, one arrives at the following conclusions. 

1. Class of equations ^ 

(a) The polynomials on the edges may be given in terms of a polynomial h{x,y;a, P), 
which is symmetric in its two first arguments. Specifically, 

. X _ f h{x,y;a,(3), \i - j\ = 1 , . 

where i ^ j and {i,j} ^ {2, 3}. 

(b) The diagonal polynomials have the same form 

hu{x,y) = h23{x,y) = G{x,y;a, (3) , (9) 
where G is symmetric in its two first arguments and in the parameters. 

2. Class of equations 

(a) The polynomials on the edges read as 

\ h2{x,y), = 2 ' ^ ^ 

where i 7^ j, {i,j} 7^ {2,3}. In the above relation, hi, h2 are biquadratic, symmetric 
polynomials of their arguments. 

^The dependence of Q on the lattice parameters a, (3 is not affected by these symmetries. That is why we omit 
them from the arguments of Q. 
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(b) The diagonal polynomials have the same form, i.e. 

hi4{x,y) = h23ix,y) = G{x,y) , (11) 

with G being symmetric. 

The most generic equation (|4]) belonging in is equation Q5, [32], for which the function Q 
has the following form. 

Q{u, X, y, z) :— aiuxyz + 02 {uyz + xyz + uxz + uxy) + 03 {ux + yz) 

+04 (uz + xy) + as [xz + uy) + 03 {u + x + y + z) + aj, (12) 

where are free parameters. 

Equation Q5 contains equations of the class ^ by choosing appropriately the free parameters 
Qi. The specific choices leading to the ABS equations HI H3 and Q1-Q4 are presented in the 
following list. 

i) Equations H1-H3 and Q1-Q3 correspond to the choice ai = 02 = 0, while the rest parameters 
are chosen according to the next table. 
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ii) Equation Q4, which is the master equation of the members of the ABS list, 5 , corresponds 
to the following choices for the parameters: 

, h n u ab{a + b) 2 /o 2 .92 w 

ai^a + b, a2 ^ -afj - ba , 03 = — 7^ + ap -{2a -)o, 

Z(a — p) 4 

ai = a(3 +ba , 05 = — - + ba -(2/3 - — )a, 

2{f3 — a) 4 

53 52 92 
ae = y ai - ^02 , 07 = — oi - 3302 , 

with = p{a), b^ — p{l3), where p{x) — Ax^ — g2X — gs- 
Remark 3.1 

Equations HI, H35=o and Ql^^g ^^'^ referred to as the discrete potential, modified and 
Schwarzian KdV equation, respectively, cf. review [ij. Equations Ql and Q35^g, are related 
to Nijhoff-Quispel-Capel (NQC) equation, [21], cf. also [3l[8]. Equation Q4 was first presented 
by Adler in [5] , as the superposition principle of the Backlund transformation for the Krichever- 
Novikov equation. 



4 Integrability aspects of the ABS equations: auto-Backlund 
transformations and Lax pairs 

An important characteristic of the ABS equations is their multidimensional consistency 3]. The 
latter means that, the equation can be consistently embedded on a three dimensional lattice, as 
illustrated in Figure [31 and, consequently, on a multidimensional lattice. 
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To be more precise, suppose that, the function u depends on a third lattice variable fc, with 
which the parameter 7 is associated. Now, each face of the cube carries a copy of the equation, 
cf. Figure [3l which involves the values of u assigned to the vertices of the face, as well as the 
corresponding parameters assigned to the edges. Thus, starting with the values of u at the black 
vertices and using the equations on the gray faces, one can uniquely evaluate the values of u 
at the gray vertices of the cube. Subsequently, the equations on the white faces provide three 
different ways to calculate the value If these different routes lead to the same result, than 

the equation is said to possess the consistency property. Moreover, if the resulting value 
is independent of U(o.o,o), then the equation fulfills the tetrahedron property. 

The ABS equations possess both of the above properties, [3. The consistency and the tetra- 
hedron properties imply that, the polynomial h related to the edges is factorized as 

h{x,y;a,f3) = k{a, f3) f{x,y,a) , 

where the function k(a,(3) is antisymmetric, i.e. 

k{(3,a) = ~k{a,P). 

Furthermore, the discriminant 

is independent of the lattice parameters. The functions /, k and G corresponding to each one of 
the equations under consideration are presented in Appendix A. 

The consistency property, which may be thought as the discrete analog of the hierarchy 
of commuting flows of integrable differential equations, [22], can be used effectively to derive 
Backlund transformations and Lax pairs, [iniinilS]- In fact, for the case of the ABS equations, 
this leads to the following result. 

Proposition 4.1 

The system 

Brf(^,u,A) ( Q("(o,o),"(i,o),«(o,o),«(i,o);a,A) = 

defines an auto-Bdcklund transformation for the ABS equation 

Q("(o,o): "(1,0)7 "(0,1), "(1,1); = 0. (14) 
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Proof 

Let us first recast equations (fT3|) and (fT4|) in terms of tlic polynomial / by using relations ([7]), 
i.e. rewrite them as 

/('«(o,o),W(i,o),a)/(u(o,o),'"(i,o),a) = /("(o,o), "(o,o), A)/(w(i,o), W(i,o), A), (15) 
/(?^(o,o),'"(o,i),/3)/(u(o,o),'"(o,i),/3) = /("(o,o),"(o,o): A)/(u(o,i),-U(o,i), A), (16) 

and 

/(M(o,o),'«(i,o),a)/("(o,i)''"(i,i)'") = /("(o,o),'"(04)'/')/("(i,o),U(i,i)'/')' (1'^) 
respectively. 

Let u be a solution of p7|) . We solve equations (fT5|) . (fT6| for the polynomials /(w(o,o)j '"(1,0)7 ck) 
and /(■U(o,o)i ^^(o,i)i /3)j respectively, and take the shift of the resulting relations in the m and n 
direction, respectively. The substitution of all these relations into (fT7|) leads to 

/('«(o,o),it(i,o)7a)/(w(o,i)>'"(i,i)7a) = /("(co), "(0,1)' /3)/("(i,o)) (18) 
which is equivalent to 

Q("(o,o)-"(i,o):'«(o,i)'"(i,i);"'/') = 0- 
Conversely, assume that u satisfies the above equation. We solve (fT5|) . (fT6|) for the polynomials 
/(m(o.o), '"(i.o) I Q;) E^nd /('ii(o,o))U(o,i),/?), respectively. The resulting expressions and their shifts, 
combined with (fT8|) . lead to equation (fT7|) . 




U2 

Figure 4: Bianchi commuting diagramm 

A direct consequence of the previous result is the next proposition. 
Proposition 4.2 (Bianchi commuting diagramm) 

Let , V? he solutions of the ABS equation {14^ , generated by means of the auto-Bdcklund 
tranformation B^; from a known solution u'^ via the Bdcklund parameters Ai and A2, respectively. 
Then, there is a third solution u^^ , which is given algebraically by 

Q{u°,u\u^,u^^;XuX2) = 0, 
and is constructed according to Figure^ 
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Proof 

Rewriting the equations of the auto-Backlund transformation in terms of the polynomial /, we 
divide the equations of u^, Ai) by the corresponding ones of (vP, v? , A2) and, the equa- 

tions of Bd A2) by the corresponding ones of B^ (-u^, m^^, Ai) . Equating the resulting 

relations, we get 

F^„(F) = G^„(G), F,yr,^{F) = G^,„(G), (19) 

where 

P = /('"('o,0)''"(0,0)'^l)/('"(0,0)''"(aO)'^l) ' 
= /('"(0,0)''"(0,0)'^2)/('"(o o):'"Jao)'^2) ■ 

From equations we get that F = G. The latter equivalently implies that 

Q{u'^,u\u\u^^;XuX2) - 0. 

Remark 4.1 

Auto-Backlund transformations different from the ones presented here, as well as some hetero- 
Backlund transformations, were given recently by Atkinson in [3]. 

The consistency property also serves in the construction of a Lax pair, through an algorithmic 
procedure [ini IB US]- Alternatively, using the fact that, a Backlund transformation may be 
regarded as a gauge transformation for the Lax pair, lO], one may construct a Lax pair for the 
ABS equations using the equations constituting B^ [33J. The two approaches are equivalent and 
lead essentially to the same result, which can be stated in the form of the following proposition. 

Proposition 4.3 

The equation of the ABS class 

Q("(o,o),'«(i,o),'«(o,i),'"(i,i);"'/?) = 

arises in the compatibility condition of the linear system 

*(i,o) = L('U(o,o),'"(i.o);a,^)*(o,o): *(o,i) = L('U(o,o), •"(0,1)! Z^' ^)*(o,o) 7 (20) 

where 

L(x\x2;a,A) = ^ , . ( '^^'^ ] , (21) 



^/c(a,A)/(a;i,a;2,a) 
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and Q — Q{x^ ^x^ ,x'^\a,\) and its derivatives (Q^i = dr^iQ) are evaluated at x^ = x'^ = 0. 



5 Infinite hierarchies of generahzed symmetries 

In this section we present the symmetry analysis of the equations belonging in the class 
while the corresponding analysis for the class was given in [30] • The conclusions of the 
symmetry analysis are summarized in the form of two propositions and, subsequently, are applied 
to equation Q5. Finally, we explicitly construct infinite hierarchies of generalized symmetries 
for all of the ABS equations, the members of which are constructed inductively using linear 
differential operators. 

Let us first present the two propositions about the symmetries of the equations belonging in 
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Proposition 5.1 

Every equation in the class admits two three point generalized symmetries with generators the 
vector fields 

^l(w(o,0),W(1.0)) 1, , 

'll,«(i.o)l'"(O,O),W(l,0)j I C«(o,o) ' 



"(1,0) ~ ""(-1,0) 2 

and 

■'"(O.O) I 



V "(0,1) - "(O.-l) 2 <°'i' V 



respectively. 
Proposition 5.2 

Let equation in the class .'P' he such that, the matrices 

, « = 2, (22) 



hi{x,y) G{x,z) G{x,w) 
hi^^{x,y) G^x{x,z) G^x{x,w) 
hi,xx{x,y) G^xx{x,z) G^xx{x,w) 



x = Q 



are invertible. Then , the generator of any five point generalized symmetry of this equation will 
be necessarily a vector field of the form 

V = a(n)v„ +6(m)v,„ + ^ij{n,m,U(^o,o))du^o,o) > 

where the functions a{n), b{m) and 'ip{n,m,U(^Q^Q^) satisfy the determining equation 

/ / N / , 1^^ J, / n2 o f ^(^(1,0), "(04)) 

(a(n) - a(n+ 1)) /ii(u(o,o), "(1,0)) Ouno) Tl T 

' \"il"(o,o), "(i,o)j 

/L( \ u , 1^^ f, / ^2a / G(m(i,o), "(o,i)) 
+ 6(m - 6 TO + 1 /i2("(o,o), "(0,1) Ou(o.i) 777 T 

V/»2("(0,0),"(0,l)) 

+G(u(i,o),U(o.i))V^(n,m,U(o,o)) + /i2("(o,o), "(o,i))'0(" + 1, "(i,o)) 
+/ii("(o,o),"(i,o))V'('^,TO+ l,"(o,i)) - + l)"(i,i)) = 0- 

Proof 

The proofs of the above propositions follow from the corresponding ones given in |30) by making 
the following changes 

h{u, X, a, P) — > hi{u, x) , h{u, x, /3, a) — > /i2(", 2:) . 

The application of the above symmetry analysis to Q5 implies that, it admits only a pair of 
three point generalized symmetries with generators the vector fields 

^ _ fel("(0,0)."(l,0)) 1 ^fel("(0,0)7"(l,0)) 

"(1,0) - "(-1,0) 2 c'u(i^o) 

^ _ fe2("(0,0),"(0,l)) 1 '9fe2("(0,0)7"(0,l)) 

"(0,1) - "(0,-1) 2 i9u(o,i) 

where 

hi{x,y) = {03 + aixy + a2{x + y)){a7 + a^xy + aeix + y)) - 
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(ae + osx + (04 + a2x)y){ae + a^x + (05 + a2x)y) , 



h2{x,y) = {a^ + aixy + a2{x + y)){aT + a^xy + aeix + y)) - 
(ae + a4X + (03 + a2x)y){a& + a^x + (04 + a2x)y) . 

Furthermore, it can be shown by a direct calculation that, the commutator of these two 
symmetry generators is a trivial symmetry, [23], i.e. the characteristic of the resulting vector field 
vanishes on solutions of equation Q5. Thus, we can write 

[v„ , v„] = . (23) 

Since Q5 contains all the ABS equations, the above commutation relation also holds for the 
corresponding symmetry generators of the latter equations. 

On the other hand, the symmetry analysis of the ABS equations, |30j, showed that, all of 
them admit a pair of three point generalized symmetries, as well as a pair of extended generalized 
symmetries, the generators of which are the vector fields 

and 

V„ - n rM 9„,„,,, - r{a) 9„ , V,„ = m - r(/3) dp , (25) 

respectively, with 

„ro] /(m(o,o),M(i,o),«) 1 , / N 

= „^ 2-'"'"(i'")'-'"(0'0)''"(l'0)'"-' ' 



i(l,0) ~ ""(-1,0) 



T,[0] /("(0.0),M(0,1)^/^) If/ m 
Rm = , 2-''"(».i)''"(0'0),U(0,l)'P)- 



^^(0,1) - ^^(0,-1) 



The form of the function r appearing in (|25p depends on the particular equation and is given in 
the following tabic. 

Equation HI 112 H3 Ql Q2 Q3 Q4 



r{x) 1 1 ~f 1 1 -f -W'^^^-92X-g3 

The generalized symmetries have been used effectually to derive reductions of partial difference 
equations to discrete analogues of the Painleve equations. The first result in this direction was 
presented by Nijhoff and Papageorgiou in [19], where the reduction of HI to a discrete analogue 
of Painleve II was given. 

The extended symmetries, which act, not only on the dependent variable u but, on the lattice 
parameters as well, are very important and can be used effectively in two different ways. The 
first way is in the construction of solutions through symmetry reductions, [331131] . Assuming that 
the function u depends continuously on the lattice parameters a, (3, one may derive solutions of 
the given difference equation which remain invariant under the action of both of the extended 
symmetries generated by the vector fields V„, Vm- In this fashion, one is led to a system 
of differential-difference equations, which can be equivalently written as an integrable system 
of differential equations. On the one hand, this differential system is related to the so called 
generating partial differential equations, introduced by Nijhoff, Hone and Joshi in [TH], cf. also 
[28l |29] . On the other hand, some of its similarity solutions reveal new connections between 
discrete and continuous Painleve equations, ^31j . 

The other way to use the extended symmetries was suggested by Rasin and Hydon in [24] . 
Specifically, it was pointed out that, the above extended generalized symmetries can be regarded 
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as master symmetries of the corresponding generalized ones. To prove this, one has to show that 
the foUowing commutation relations hold 



V v["l 



^ 0, 



V v["l 



[0] 



= 0, 



as well as similar relations for the generators V^, . It should be noted that, the commutators 
lead to trivial symmetries because of the relation (^5]) . e.g. [V„,Vm] = 

Now, writing out explicitly the commutators [V„,v.[?'] and [V„i,Vm], one arrives at 



V v[°ll - Rl^la 



V 



(26) 
(27) 



where 



(^„R[°1) - (^i-^'R™) - r(a)5„) R™ 

_ /(^(0,0). ""(1,0); Q)/(m(0,0): «)(m(2,0) - '"(-2,0)) 

(W(1,0) - U(-1,0))^(W(2,0) - '"(0,0))("(-2,0) - "(0,0)) 

rW = ((^,„rm) - (^ir'^Rj^i) - r(/?)a;3) r™ 

_ ./(w(0,0),'"(0,l)'/')/('"(0,0),'«(0,-l),/3)('"(0,2) -W(0-2)) 



(28) 



(W(0,1) - M(0,-l))^("(0,2) - W(0,0))(W(0 -2) - U(0,0)) 

Using definitions (|26ll27p and the Jacobi identity, it can be verified straightforwardly that: 



(29) 



n ' m 



ni ' n 



v[il V 



= 0. 



On the other hand, the commutators 



* n ' * n. 



771 ' 



= 



can be verified by using pi]) . (|26ll29p and the properties of polynomial /, i.e. it is biquadratic 
and symmetric. 

The above analysis implies that, the symmetry generators V„, are master symmetries of 
and Vm , respectively. Consequently, an infinite hierarchy of generalized symmetries can be 
constructed in this fashion, the members of which are defined inductively: 



"(0,0) 



[k] 



, fc = 0, 1,... and i = n, m . 



(30) 



The characteristic R^'"' involves the values of u at {2k + 3) points in the i direction of the lattice 
and is determined by applying successively the linear differential operators 



oo 

= ^ £ (^i^'R[°l)a.,,,„, - r(a)9„, 

oo 



(31) 
(32) 
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on ' , I.e. 



Rf ' = ^.''Rf , fc = 0, 1 . . . and i ^ n,m. (33) 



Remark 5.1 

The linear operators rnay be regarded as recursion operators for equations Q3 and Q4. 

This foUows from the fact, [30] . that, these equations admit only the symmetries generated by 
the vector fields given in ([24|) . 

Remark 5.2 

It is worth mentioning some previous results on higher order generalized symmetries for the ABS 
equations. Specifically, a procedure for the construction of hierarchies of generalized symmetries 
was presented in [13], where the relation of the ABS equations to Yamilov's discrete Krichever- 
Novikov equation (YdKN), [34l[35], was explored. On the other hand, particular results about 
HI and Ql^^o were presented in [11] and [12], respectively, the derivation of which based on the 
associated spectral problem. The advantage of our approach compared to the above outcomes is 
that, it leads straightforwardly to explicit expressions for the higher order generalized symmetries 
for all of the ABS equations. 

The correspondence of the generators vp' to group of transformations, [23], leads to hierarchies 
of integrable differential-difference equations, the first members of which are special cases of YdKN 
equation, [13], cf. also [5 . 

Dropping the one of the two indices corresponding to shifts with respect to n or m, and 
denoting by a the corresponding lattice parameter, these hierarchies have the form 

du ( f(u,ui,a) 1 „ , , „ , 

dcfc V "1 - ""-1 2 J 

where is the corresponding group parameter, 

^= E ^ (^^'^ ( ^I'^-u"! - - ^(^)9a (35) 

and is the shift operator in the corresponding direction. 
The corresponding ABS equation, written now in the form 

Q{u,ui,u,ui; a, X) = 0, 



is an auto-Backlund transformation of equations p4p . [13 , while the latter equations admit the 
following Lax pair 

d^" 

*i = L(w,ui;a,A)*, — = Nfc«', A: = 0, 1, . . . . (36) 

In the above relations, the matrix L(u, ui; a, A) is given by (j2ip and 

1 1 
X - - 

ui — U-i 2 



Nfc ( X - i X^„, ) , fc = 0, 1, . . . , (37) 



where 



X :— — f{u,ui,a) (L(u, wi; a. A)) ^ (?„L(u, wi; a. A) . (38) 
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6 Conclusions and prespectives 



We have presented some recent results about difference equations, their integrabihty aspects and 
their underlying symmetry structure. We have used as a representative example the equations 
of the ABS classification |3]. The multidimensional consistency property possessed by the ABS 
equations provides the means to investigate in a systematic way Backlund transformations and 
Lax pairs for all of them. In the same fashion, integrabihty aspects of consistent systems of 
difference equations have been derived, e.g. the discrete Boussinesq system, [26l [27], and the 
elliptic generalization of the discrete KdV equation, [20] . 

The extended generalized symmetries admitted by the ABS equations are proven to provide us 
an effective tool. They play the key role of master symmetries leading to the explicit construction 
of infinite hierarchies of generalized symmetries. Moreover, the same symmetries can be used 
for the construction of similarity solutions, which were called continuously invariant solutions 
in [31', revealing a new link among difference and differential equations. Moreover, continuously 
invariant solutions of special type have exposed a novel connection among discrete and continuous 
Painleve equations. It would be interesting to study other systems of difference equations from 
the point of view of extended generalized symmetries and continuously invariant solutions, such 
as the discrete Boussinesq system, for which their connection to integrable differential equations 
was established by different means in [26l [27] . 

A The characteristic polynomials of the ABS equations 

HI /(m, x,a) = 1 

G{x,y) = (x-y)^ 
k{a, (3) = P — a 

H2 /(m, X, a) = 2{u + x + a) 

G{x,y) = {x~yY-{a-Pf 
k{a, (3) = 13 — a 

H3 /(m, X, a) = ux + a5 

G{x, y) = {ya - x(3){yl3 - xa) 

Ql f{u,x,a) — {{u — x)"^ — 5^) / a 
G{x,y)^a(3{ix-y)^-ia-P)H') 
k{a,P) = -a(3{a-P) 

Q2 /(u, X, a) = {{u - xf - 2a'^(u + x) + a^)/a 

G{x, y) = aP {{x - yf ~ 2{a - pY{x + y) + (a - pf) 
k{a,P) = -aP{a-P) 

Q3 /(u, X, a) = 4^(~Li) i-iaix - au){ax - u) + {a^ - ifS"^) 

Gix,y) = ^"'-'^^f-'^ (4a/3(«y - Px)iPy - ax) + {o? - P^)5^) 
k{a, P) ^ (a2 _ ^2)(^2 _ i)(^2 _ 1) 

Q4 /(u, X, a) = {(ux + a{u + x) + 52/4)'^ — {u + x + a){4:aux — g^)) /a 

G(x, y) = {aixy + a2{x + y) + ai){aixy + aQ{x + y) + 07) 

- {a2xy + 03?/ + a^x + a(,){a2xy + 032: + a^y + a^) 
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^\ _ ab(a^6+aft^ + [l2a/3'^-g2(a+2/3)-3g3]a+[l2/3a''-g2(/3+2a)-3fl3]b) 
K(a,P) - 4(a-/3) 
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